Confirmation of quark-hadron duality in the neutron F2 structure function 
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We apply a recently developed technique to extract for the first time the neutron F% structure 
function from inclusive proton and deuteron data in the nucleon resonance region, and test the 
validity of quark-hadron duality in the neutron. We establish the accuracy of duality in the low- 
lying neutron resonance regions over a range of Q 2 , and compare with the corresponding results 
on the proton and with theoretical expectations. The confirmation of duality in both the neutron 
and proton opens the possibility of using resonance region data to constrain parton distributions at 
large x. 



The quest to understand the strong interactions at in- 
termediate energies, and particularly the transition from 
quark-gluon to hadron degrees of freedom, is one of the 
main outstanding challenges in modern nuclear physics. 
Considerable attention has been focused recently on the 
"duality" between quark and hadron descriptions of ob- 
servables in electron-hadron scattering. A classic exam- 
ple is the finding [l[ that inclusive structure functions in 
the region dominated by the nucleon resonances on av- 
erage resemble the structure functions measured in the 
deep inelastic scattering (DIS) region at higher energies. 

With the availability of high-precision data from Jef- 
ferson Lab and elsewhere, duality has now been firmly es- 
tablished for the proton F 2 and Fr. structure functions 0- 
H[, and exploratory studies in spin-dependent and semi- 
inclusive scattering have provided tantalizing glimpses of 
the flavor and spin dependence of duality (for a review 
see Ref . @ ) . A complete picture of the workings of du- 
ality in the nucleon can only be constructed, however, 
with information on duality in the neutron, on which lit- 
tle empirical data exists. 

Calculations based on quark models point to intrigu- 
ing differences between duality in the proton and neutron 
@ , and some arguments even suggest that duality in the 
proton may be due to accidental cancellations between 
quark charges, which do not occur for the neutron 0. 
Confirmation of duality in the neutron would therefore 
firmly establish that the phenomenon is not accidental, 
but rather a robust feature of nucleon structure func- 
tions. More generally, understanding the transition be- 
tween the resonance and DIS regions can lead to better 
constraints on parton distribution functions (PDFs) at 
large momentum fractions x, by allowing data at lower 
final state hadron masses W to be used in global PDF 
fits [I,©]. Precise knowledge of large- a; PDFs, which are 
currently poorly constrained, is vital in searches for new 
physics beyond the Standard Model [lfj, for instance, as 



well as in neutrino oscillation experiments [11 



In this Letter we use a recently introduced technique 
[l2T ] to extract for the first time the neutron F 2 structure 
function from proton (p) and deuteron (d) F 2 data in the 
resonance region over a range of photon virtualities from 
Q 2 = 0.6 to 6.4 GcV 2 . The new method is based on 



an iterative approach in which the nuclear corrections 
are applied additively, and has been found to accurately 
reproduce neutron structure functions of almost arbitrary 
shape in both the DIS and resonance regions [l3 |. 

The extraction of reliable neutron information from 
deuterium data requires a careful treatment of nuclear 
effects [l3|], and we use the latest theoretical develop- 
ments which allow the deuteron structure function to be 
analyzed in both the resonance and DIS regions, at both 
low and high Q 2 . In the weak binding approximation 
the deuteron if structure function can be written as a 
sum of smeared proton and neutron structure functions 
F 2 {N = p, n), and an additive term which accounts for 
possible modification of the structure functions off-shell 



if = if + if + 5 {oS) F* . 



(1) 



The smeared nucleon structure functions are given by 
convolutions of the nucleon light-cone momentum dis- 
tribution in the deuteron, fiy/d, and the bound nucleon 
structure functions [12, E3] , 
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where the symbol C§> denotes a convolution. The nucleon 
momentum distribution (or smearing) function /jy/d ac- 
counts for the effects of the nucleon's Fermi motion and 
binding, including finite-Q 2 corrections [H, E3, and is 
taken to be identical for the proton and neutron. The off- 
shell correction S^° S ^F 2 has been found in several models 
to be typically of the order 1-2% for x < 0.9. 

To account for the quasi-elastic (QE) tail in the 
deuteron data the clastic nucleon contribution is smeared 
using the same /jy/d- Subtracting from the deuteron if 
the QE contribution, together with the off-shell correc- 
tion and the smeared proton if, one obtains an effective 
smeared neutron structure function F 2 and then solves 
Eq. ([2]) for the neutron. 

The nuclear effects are parametrized by an additive 
correction 1121. 



NFS + 5f® if . 



(3) 
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0, [ljl and half the difference in the results assigned 
as a systematic uncertainty. 

The stability of the iteration method relies on the avail- 
ability of relatively smooth data, especially for deuterium 
(irregularities in proton data are smoothed out by the 
smearing). This is critical at large x where the structure 
functions are small, and discontinuities could even render 
the extracted neutron results negative. It is particularly 
important that the QE contribution to the deuteron F$ 
be accounted for in the analysis, and we model this using 
the same smearing function, /jv/di an d nucleon form fac- 
tors from Refs. [HI, This is found to provide a good 
description of the QE peak as a function of Q 2 . 

An example of the extracted neutron F£ structure 
function is displayed in Fig. [T] for Q 2 = 1.7 and 5 GeV 2 , 
together with the input proton and deuteron data (the 
complete data set will be shown in Ref. [HI). The start- 



ing value of the neutron for the iteration was F : 



n(0) 



F p 
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FIG. 1: Extracted neutron F% structure function at Q 2 = 1.7 
and 5 GeV 2 , together with proton and deuteron data, and 
the reconstructed deuteron (total uncertainties are system- 
atic and statistical errors added in quadrature). The proton 
and neutron data are compared with the global QCD fit from 
Alekhin et al. ||. The dependence of the iteration on the 
initial value is illustrated in the insert (see text). 



and the deuteron F% reconstructed from the proton and 
extracted neutron was found to be in good agreement 
with the data after two iterations. The spectrum of the 
F£ structure function in the resonance region displays 
similar characteristics as observed from the proton spec- 
trum: one finds three resonant enhancements which fall 
with Q 2 at a similar rate as for the proton. 

To check that the extracted neutron structure function 
does not depend on the starting value of the iteration, the 
extraction procedure was repeated assuming a different 
boundary condition, F^ ^ = Ff/2. The difference be- 



where Sf gives the finite width of the smearing function tween the two results A = [F£(F< 
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and Af its normalization. The F^ structure function is 
then extracted using an iterative procedure [12], which 
after one iteration gives 
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starting from a first estimate F^^, and iterated until 
convergence is reached. The robustness of this method 
and its ability to reliably estimate errors on the ex- 
tracted neutron function were investigated extensively for 
smooth functions in Ref. 12 |. Since the smearing func- 



tion is sharply peaked, the convergence of this method 
is typically extremely fast, requiring only one or two it- 
erations before the F$ function reconstructed from the 
extracted F£ matches the original data to within exper- 
imental uncertainties. 

In this analysis we use proton and deuteron Fi data 
from JLab experiment E00-116 Q] and SLAC experi- 
ments E49a6 and E49al0 [l7| . The former span the high- 
Q 2 region, 4.5 < Q 2 < 6.4 GeV 2 , while the latter cover 
the lower Q 2 range, 0.6 < Q 2 < 2.4 GeV 2 , providing a 
total of 514 data points. Because the extraction method 
requires proton and deuteron F 2 data at fixed Q 2 , the 
centering of the data at the same Q 2 was made at the 
cross section level using the p and d fits from Ref. [l8j |. 
To test the sensitivity of the results to the choice of the 
bin-centering fit, an additional fit was used for each target 



F p /2)]/a(F2), normalized by the total F£ uncertainty 
c(i<2 i ), is shown in the insert of Fig. [1] after two itera- 
tions. One finds an almost Gaussian distribution cen- 
tered around (the mean of the distribution is around 
—0.07) with a width well within the typical total uncer- 
tainty of F£ . In fact only 6% of the total number of data 
points lie outside of a 2a range. More extreme boundary 
conditions, such as F^^ = 0, do not alter the character- 
istics of the extracted F^ structure function spectrum, 
with the resonant structures already visible after just 1 
iteration. On the other hand, as discussed in Ref. fl2| . 
more iterations are needed for poor choices of initial val- 
ues, which increases the scatter of data points if the deu- 
terium data in particular display any nonuniformitics. 

The effect of the off-shell correction was taken 

into account using the model of Ref. [l6j|, which gives 
~ —1.5% correction over most of the x range considered, 
and was argued to provide an upper limit on the correc- 
tion. The F^ data are corrected by subtracting half of 
the off-shell correction from [l6[ and assigning a 100% 
uncertainty. When propagated into the F£ uncertainty 
this was found to contribute less than 2% to the total 
error. 

In Fig. Q] we also show F% and F% from global QCD 
fits to DIS data (with W 2 > 4 GeV 2 ) from Alekhin et al. 
Q, which illustrates the striking similarity between the 
QCD fit and the resonance data, reminiscent of Bloom- 
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FIG. 2: Truncated neutron moments AfJ ("data") in various 
resonance regions (1 st , 2 nd , 3 rd and W < 2 GeV) relative to 
the QCD fit from Alekhin et al. ("theory"). 

Gilman duality To quantify this duality we consider 
ratios of "truncated" moments Mi [Hf , 

A/ 2 (Q 2 ,A.t) = / dxF 2 (x,Q 2 ), (5) 

J Ax 

in the resonance region for specific intervals Ax. Follow- 
ing previous proton data analyses @, Q , we consider the 
regions 

• 1 st resonance region — > W 2 G [1.3, 1.9] GeV 2 

• 2 nd resonance region — > W 2 G [1.9, 2.5] GeV 2 

• 3 rd resonance region — > W 2 G [2.5, 3.1] GeV 2 

as well as the entire resonance region 1.3 < W 2 < 
4 GeV 2 . At a given Q 2 , the lowest- W (A resonance) re- 
gion corresponds to the highest- a; range, and for a fixed 
W interval the larger the Q 2 , the higher the x. 

The ratio of the truncated moments of the resonance 
data to the global QCD fit , computed over the same x 
range, is shown in Fig.[5]as a function of Q 2 . Globally, the 
agreement between the QCD fit and the resonance data 
is quite remarkable, with deviations of < 10% observed 
over the entire Q 2 range. Locally, in the individual res- 
onance regions the deviations are generally < 15 — 20%, 
somewhat larger only in the 1 st resonance region at the 
largest Q 2 . This is not surprising given the fact that the 
A region at Q 2 = 6.4 GeV 2 covers the highest-x regime 
studied, x ~ 0.9, where the QCD fit is mostly beyond its 
limit of applicability. 



FIG. 3: Ratio of truncated neutron to proton moments 
M2/M2 in various resonance regions as a function of Q 2 , 
compared with global QCD fits from Alekhin et al. [s| and 
MSTW [13. 



The isospin dependence of duality can be studied by 
comparing the truncated neutron moments with the anal- 
ogous proton moments. The ratio of these is displayed in 
Fig. [3] as a function of Q 2 for the various resonance re- 
gions, and compared with global QCD fits from Alekhin 
et al. M a nd from MSTW [25| , corrected for target mass 
effects [26|. The MSTW fits are shown for Q 2 > 2 GeV 2 , 
which corresponds to their approximate limit of validity. 

The ratios show good agreement with the data, with 
the exception of the A region which is somewhat under- 
estimated. Since the proton and neutron transitions to 
the A are isovector, the resonant contributions should be 
identical; on the other hand, the DIS structure functions 
in the A region are expected to be rather different, with 
F-2 <C F%, so that violation of duality here is expected 
to be strongest. In addition, the QCD fits are least con- 
strained in this region due to the scarcity of large- a; DIS 
data. This is especially the case for the MSTW fit [H| 
which limits the data sets to W 2 > 15 GeV 2 . 

The M2/M2 ratios at fixed Q 2 are shown in Fig. @] as 
a function of x for the three resonance regions, compared 
with the QCD fits as in Fig. [3] The global fits offer 
a good description of the 2 nd and 3 rd resonance region 
data, revealing clear evidence of duality down to Q 2 as 
low as 0.6 GeV 2 . The fits underestimate the A-region 
ratios and this trend becomes more pronounced as one 
moves to larger Q 2 (> 4 GeV 2 ) and larger x. The Alekhin 
et al. fit @ offers a better description at large x, which 
is likely due to its inclusion of lower- IF, lower-Q 2 data. 

Our results can be compared with quark model ex- 
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FIG. 4: As in Fig. but as a function of £ for fixed Q 2 



pcctations for the isospin dependence of duality, which 
predict systematic deviations of resonance data from lo- 
cal duality. Assuming dominance of magnetic coupling, 
the proton data are expected to overestimate the DIS 
function in the 2 nd and 3 rd resonance regions due to the 
relative strengths of couplings to odd-parity resonances, 
especially those in the quark spin-i octet Q. The neu- 
tron data are predicted to lie below the DIS curve in 
the 2 nd resonance region due to the small coupling to 
octet states with spin i, but have larger couplings to 
odd-parity quark spin-| octet states. Remarkably, the 
neutron data do indeed underestimate the global F"£ fits 



in the 2 nd resonance region, just as theproton data were 
found to exceed the global F^ fits 0, 0]. Moreover, the 
similarity between the truncated M£ moments in the 
W 2 < 4 GeV 2 and DIS regions strongly suggests that the 
resonance cancellations in the proton are not accidental 
0] , but rather form a systematic pattern which dramat- 
ically reveals itself through the Bloom-Gilman duality 
phenomenon. 

In conclusion, we have extracted the neutron structure 
function for the first time in the resonance region 
from inclusive proton and deuteron data. Our compar- 
isons of empirical truncated moments to those extracted 
from global QCD fits to high-W 2 data show clear signa- 
tures of Bloom-Gilman duality, with better than 15 — 20% 
agreement in the 2 nd and 3 rd resonance regions, and 
less than 10% deviations when integrated over the en- 
tire W 2 < 4 GeV 2 region. The confirmation of duality 
in the neutron establishes that the phenomenon is not 
accidental, but is a general property of nucleon structure 
functions. Our findings suggest that averaged resonance 
data could in future be used to constrain the large- a; be- 
havior of global QCD fits by relaxing the W 2 cuts on 
data down to the 2 nd resonance region, W 2 ~ 1.9 GeV 2 . 
This could also have significant impact for searches for 
new physics beyond the Standard Model at colliders and 
neutrino oscillations experiments. 
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